Collective dynamics and expansion of a Bose-Einstein condensate 

in a random potential 
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We investigate the dynamics of a Bose-Einstein condensate in the presence of a random potential 
created by optical speckles. We first consider the effect of a weak disorder on the dipole and 
quadrupole collective oscillations, finding uncorrelated frequency shifts of the two modes with respect 
to the pure harmonic case. This behaviour, predicted by a sum rules approach, is confirmed by the 
numerical solution of the Gross-Pitaevskii equation. Then we analyze the role of disorder on the one- 
dimensional expansion in an optical guide, discussing possible localization effects. Our theoretical 
analysis provides a useful insight into the recent experiments performed at LENS [J. Lye et al, 
Phys. Rev. Lett. 95, 07040f (2005); C. Fort et aZ., |cond-m at/0507144|. 



I. INTRODUCTION 

The investigation of Bose-Einstein condensates (BECs) 
in the presence of disorder is rapidly becoming a central 
topic in ultracold atom physics (HHHSlSlESIi- 
Bosonic systems in disordered potentials have been ex- 
tensively investigated in the recent past, both experimen- 
tally and theoretically Q. Experiments with superfluid 
He in porous materials have demonstrated the suppres- 
sion of superfluid transport and the critical behavior at 
the phase transition in presence of disorder |Toj . From 
the theoretical point of view a rich variety of phenom- 
ena is expected to occur in these systems, among which 
the most fascinating are Anderson localization, initially 
proposed in the context of electron transport in disor- 
dered solids and later predicted and observed for 
non- interacting wave phenomena such as light 0, ^| , 
and the quantum transition to the Bose glass phase that 
originates from the interplay of interactions and disorder 



The demonstrated capability of using BECs as ver- 
satile tools to revisit condensed matter physics 0, as 
for example the transition from superfluid to Mott in- 
sulator jig, suggests that these are also promising tools 
to engineer disordered quantum systems Q, 0, U ■ Re- 
cently, effects of disorder created by a laser speckle have 
been observed on the dynamics of a BEC, including un- 
correlated shifts of the quadrupole and dipole modes 3 
and localization phenomena during the expansion in a 
one-dimensional (ID) waveguide UE|. Effects of disor- 
der have also been observed for BECs in microtraps as a 
consequence of intrinsic defects in the fabrication of the 
microchip |l7l.lTs|. 



In Refs. [J, |6| we have shown that the main features 
observed in that experiments can be explained within 
the Gross-Pitaevskii (CP) theory. In this paper we re- 
port a detailed analysis and discussion of the theoretical 
approach used, comparing the effects of different kinds 
of random potentials. We also make a systematic com- 
parison with the case of a periodic lattice with spacing 



of the order of the length scale of disorder. This helps to 
discriminate the effects due to the particular realization 
of the random potential from those that are intimately 
connected to the disorder. 

We show that in the presence of a weak disorder the 
dipole and quadrupole modes of a harmonically trapped 
condensate are undamped in the small amplitude regime, 
whereas a superfluid breakdown may occur for larger os- 
cillations. In the first case the two modes are character- 
ized by uncorrelated frequency shifts, both in sign and 
amplitude, that depend on the particular realization of 
the perturbing potential. The average features however 
do not depend crucially on the particular kind of disorder, 
but still evidence significant differences with the periodic 
case. We also show that the localization effects observed 
during the expansion in a ID waveguide are mainly due 
to a classical trapping into single wells or between barri- 
ers of the random potential. The qualitative behavior in 
this case is very similar to that of a periodic system. 

The paper is organized as follows: we start in Sect, [B] 
by describing the system and the various kinds of disor- 
der considered. Then in Sect. IIIII we discuss the effect 
of the random potential on the dipole and quadrupole 
collective oscillations of the system by means of a sum 
rules approach and the direct solution of the GP equa- 
tion. Next, in Sect. II VI we address the role of disorder on 
the BEC expansion in a ID waveguide by analyzing the 
results of the GP calculations in terms of the quantum be- 
haviour of a single defect (well/barrier) of the potential. 
A detailed description on the numerical characterization 
of the random potentials is reported in the Appendix. 



II. DESCRIPTION OF THE SYSTEM 

In this paper we will consider the case of an elon- 
gated condensate confined in a cylindrically symmetric 
harmonic potential 
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FIG. 1: (Color online) Typical shape of the potentials consid- 
ered in this paper: (a) (blue detuned) speckles (red speckles 
would have the same shape but reversed sign), (b) gaussian 
random, (c) periodic. In all the cases the correlation length 
is l c — 10/xm. 



and subjected to an additional random potential Vr{z) 
along the axial direction. The latter is characterized by 
the correlation length l c and the amplitude Vo, and can 
be written as Vr{z) = Vqv(z) with the distribution of 
intensities of v(z) being normalized to unit standard de- 
viation. Here we will address three kinds of disorder: two 
corresponding to a laser speckle potential ±Vs(z), where 
the ± indicates whether it is red or blue detuned (that is, 
the potential can be attractive or repulsive) , and another 
generated by a gaussian random potential Va(z). A de- 
tailed description on how the potentials are constructed 
and characterized is reported in Appendix lAl 

In some cases it will also be useful to compare the ef- 
fect of disorder to the case of a periodic lattice Vl{z) — 
VqVl(z). For this purpose a suitable choice is a sinu- 
soidal potential vl{z) = 2 \/2 sin 2 (nz/ 21 c ) with intensity 
normalized as before and whose wavevector is chosen to 
match the correlation length of the random potential, 
as discussed in Appendix ^] The typical shape of the 
various potentials for a correlation length l c = 10/xm is 
depicted in Fig^ 



III. COLLECTIVE EXCITATIONS 

Let us start by discussing the effect of a random po- 
tential on the collective excitations of the system, con- 
sidering in particular the dipole and quadrupole modes. 
First we will consider the regime of weak disorder and 
small amplitude oscillations, by comparing the predic- 
tion of sum rules with the numerical solution of the GP 
equation, analyzing in more details the theoretical de- 
scription of the experiments reported in Q . At the end 
of the section we will also briefly discuss the possibility 
of a superfluidity breakdown that may occur for larger 
amplitude oscillations. 



A. Sum rules approach 

A powerful tool for characterizing the collective fre- 
quencies of the system is the sum rules approach 
Within this approach an upper bound for the frequencies 
of the low lying collective excitations of a many-body sys- 
tem is given by 



1 m 3 
ft 2 mi 



(2) 



where the moments m t arc defined via the following com- 
mutators 



mi = {[F,[H,F]]) 

m 3 = ([[F,H], [H, [H,F]]]) 



(3) 
(4) 



between the many-body hamiltonian H and a suitable 
excitation operator F that is chosen as follows 
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a being a variational parameter. In our case the hamil- 
tonian can be written as 
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where the interaction strength g is related to the inter- 
atomic scattering length a by g — 4i:h 2 a/m, m being the 
atomic mass. 

In case of the harmonic potential Vho alone (unper- 
turbed case) the dipole and quadrupole collective fre- 
quencies have the well known expressions lud — to z for 
the dipole mode along z, and ujq — yjb/2uj z for the 
quadrupole mode in case of an elongated condensate in 
the large N Thomas-Fermi (TF) limit. 

Let us now discuss the effect of a shallow random po- 
tential Vr{z). Treating the Vr as a small perturbation, 
and writing u> 2 = lo 2 + 5 we get 
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where the averages (")o are calculated on the unper- 
turbed ground state, and the second line is obtained as- 
suming a strongly elongated condensate. 

The above equations imply that in general the shifts 
of two frequencies are uncorrelated and depend on the 
particular shape of the perturbing potential and on its 
relative position respect to the harmonic potential. 

To show how this works in a particular example we 
consider here the typical parameters of the LENS exper- 
iment in Q : frequencies oj z = 2ir x 9 Hz and oj± = 2n x 90 
Hz, total number of atoms N = 1 x 10 5 , Vo = 2.5 Tilo z 
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FIG. 2: (Color online) (a) Column density of a typical ground 
state configuration in the presence of a red-detuned speckle 
potential (Vb = 2.5hu> z ). (b) Solid line: profile of the com- 
bined potential Vr + Vho (dashed line). 



and l c = 10/im. A picture of the total potential resulting 
from these parameters and of the corresponding ground 
state is shown in Fig. [5] 

In Fig. |31 we show the frequency shifts Auj = u — w — 
S/2ujq and their statistical distributions respectively for 
100 and 1000 different realizations of the speckle poten- 
tial. The picture shows that the dipole and quadrupole 
shifts are uncorrelated, in contrast to what happens in 
case of a pure harmonic potential or in the presence of 
a periodic potential |2l|. This behavior does not depend 
on whether the speckles are red or blue detuned (accord- 
ing to Eqs. |H1 and this corresponds to a change of sign, 
and therefore the statistical properties in Fig. remain 
unchanged). We have also verified that the behavior is 
essentially the same also in case of a gaussian disorder, 
see Fig. 0] 

It is useful also to comment on the behaviour in the 
presence of a periodic potential. When the wavelength 
n/q = 21 c of the potential is much smaller that the axial 
extent of the condensate one can apply the Bloch picture 
and resort to the effective mass approximation. As stated 
above, this yields the same renormalization for both the 
dipole and quadrupole freq uencies u) — ^m/m*^, m* 
being the effective mass [21| . Differently, in the case con- 
sidered here (l c — 10/im) the condensate extends over 
only few wells of the periodic potential and the Bloch 
picture cannot be applied. In this case the sum rules 
approach predicts a sign correlated shift for the two fre- 
quencies, whose magnitude however still depends on the 
relative position between the condensate and the periodic 
potential. Therefore, although Aluq/ Awrj ~ constant, 
the fact that the ratio loq/lod depends at first order on 
the difference Awq — Aw^ eventually yields an uncorre- 
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FIG. 3: (Color online) Top: dipole (filled circles) and 
quadrupole (empty circles) frequency shifts for 100 different 
realizations of the speckle potential as obtained from the sum 
rules. Bottom: their probability distribution P for 1000 real- 
izations (left:dipole, right:quadmpole); the dashed lines are a 
gaussian fit. 
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FIG. 4: (Color online) Top: dipole (filled circles) and 
quadrupole (empty circles) frequency shifts for 100 different 
realizations of the gaussian random potential as obtained from 
the sum rules. Bottom: their probability distribution P for 
1000 realizations (left:dipole, right: quadrupole); the dashed 
lines are a gaussian fit. 



lated renormalization of the two frequencies (see Fig.[3J). 



B. GP dynamics 

The prediction of the sum rules can be directly com- 
pared with the solution of the GP equation 



ihdttp = 



2m 



4> 



(10) 



by exciting the collective modes with a sudden displace- 
ment of the harmonic trap (for the dipole) or a change 
of the axial trapping frequency (for the quadrupole). 
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FIG. 5: (Color online) Ratios between quadrupole and dipole 
shifts (top) and the corresponding frequencies (bottom) for 
speckle (filled circles) and periodic (empty circles) potentials 
with l c = lOpim. Note that although the shifts for the periodic 
case are correlated in sign the two frequencies are uncorrelated 
and depend on the actual position of the condensate in the 
periodic potential. 
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FIG. 6: (Color online) Dipole (top) and quadrupole (bottom) 
frequency shifts as obtained from the GP equation (squares) 
and compared with the sum rules predictions (circles), for ten 
different speckle realizations. The agreement is remarkable 
0. 



The results from some sample realization of the ran- 
dom potential are shown in Fig. [fj] (they correspond to 
the first ten realizations in Fig. [^J. The dipole oscil- 
lations are induced after a displacement Az = 5 /im of 
the harmonic potential, corresponding to an oscillation 
of the order of 10% of the axial size of the condensate. 
For the quadrupole mode, an oscillation of the same or- 
der of magnitude is obtained by releasing the condensate 
from a tighter trap of axial frequency uj' z — 1.1cj z . 

In this regime of small amplitude oscillations the so- 
lution of the GP equation shows that the condensate 
oscillates coherently with no appreciable damping on a 
timescale of several oscillations. The corresponding fre- 
quencies show a remarkable agreement with the sum rules 
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FIG. 7: (Color online) Dipole oscillations in the speckle po- 
tential (Vo = 2.5hu) z ) for three different displacements Az of 
the harmonic potential (Az = 6, 15, 24 jj,m). The insets show 
a density plot for Az = 24/im at t — 75 ms (left), Az — 15/xm 
at t — 250 ms (center), and Az — 6pm at t = 500 ms (right). 



predictions, regarding both the sign and order of magni- 
tude of the shift, as shown in Fig. El 23]. As mentioned 
above, these features have been observed in the experi- 
ment reported in Q. 

We have also explored the behavior of the system for 
larger amplitudes in case of dipole oscillations, as shown 
in Fig. [7| As the amplitude is increased the frequency 
shift reduces owing to the fact that the condensate expe- 
riences outer regions of the harmonic potential where the 
effect of the random potential is negligible, yielding an 
average frequency that is closer to the unperturbed value. 
However, as the center of mass velocity increases, the os- 
cillations also gets damped due to the presence of the 
speckle potential that acts as an external perturbation 
or roughness of the medium. In this regime the conden- 
sate develops short wavelength density modulation that 
may eventually lead to a breakdown of the superfluid 
flow, as shown in the left and center insets in Fig. [7\ (re- 
spectively for Az — 24/im at t = 75 ms and Az — 15/im 
at t — 250 ms). The rightmost inset demonstrates in- 
stead that for small displacements the condensate re- 
mains coherent even after several oscillations (Az = 6/zm 
at t = 500 ms). 



IV. EXPANSION IN A WAVEGUIDE 

Let us now consider the expansion of the condensate in 
an optical waveguide, in the presence of disorder. In this 
case we will refer to a second experiment performed at 
LENS 6] . Similar experiments have also been performed 
by D. Clement et al. |{|. The condensate is initially 
confined in an optical harmonic trap of frequencies to z — 
2ir x 30 Hz and ojj_ = 2ir x 300 Hz in the presence of 
a speckle potential of intensity Vo = 0.2 [Itf, {^tf — 
87hu> z is TF chemical potential of the condensate in the 
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FIG. 8: (Color online) (a)-(d) Density profiles of the con- 
densate (red continuous line) during the expansion in the 
waveguide in the presence of a red-detuned speckle poten- 
tial (shown in frame (e)), for different times (t = 0, 25, 50, 75 
ms, from (a) to (d)), compared with the free expansion case 
(blue dashed line). Notice that the y-axis scale changes from 
(a) to (d). 
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FIG. 9: (Color online) (a)-(d) Density profiles of the con- 
densate (red continuous line) during the expansion in the 
waveguide in the presence of a blue-detuned speckle poten- 
tial (shown in frame (e)), for different times (t = 0, 25, 50, 75 
ms, from (a) to (d)), compared with the free expansion case 
(blue dashed line). Notice that the y-axis scale changes from 
(a) to (d). 



optical harmonic trap). The condensate is prepared in 
the ground state of the combined potential, and then let 
expand through the waveguide by switching off the axial 
trapping. 

Owing to the strong radial confinement the expansion 
of the system can be conveniently described by the non- 
polynomial Schrodinger equation (NPSE) pjj . that can 
be written in a compact form as 



h 2 „ 
Am 
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(11) 

with cr 2 = y/1 + 2aN\ip\ 2 . 

In Figs. 1811 II we show the density profiles of the con- 
densate at different times during the expansion in the 
waveguide, for different choices of the random potential. 
For comparison we also show the corresponding profiles 
in case of a free expansion in the waveguide. 

Let us discuss the figures by starting from the red- 
detuned speckles (the case of Ref. 6]) in Fig. |SJ In 



this case the dynamics is characterized by an almost free 
expansion of the lateral wings of the condensate, whereas 
the central part remains localized in the deepest wells of 
the potential. This behavior can be easily explained by 
recalling that in the TF regime and in the absence of 
disorder the velocity field has a linear dependence on z, 
v(z,t) = z\ z (t)/\ z (t) (X z (t) is a scaling parameter |2^). 
indicating that the most energetic atoms reside at the 
edges of the condensate whereas the atoms close to the 
center have a nearly vanishing velocity. The presence of a 
weak disorder does not modify substantially this picture. 
This means that the outer part of the condensate can 
be sufficiently energetic to pass over the defects of the 
potential expanding as in the unperturbed case, whereas 
the central part remains partially localized in the initially 
occupied wells (see the two density peaks in the center of 
the figures). A closer look to Fig.[5]shows also that in the 
intermediate region the density distribution shows peaks 
that are instead in correspondence of the maxima of the 
potential as a consequence of the acceleration acquired 
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FIG. 10: (Color online) Density plots of the function q(k, a, 0) 
(see text) for a potential well ((a)-(b)) and a barrier ((c)-(d)), 
as a function of the incident momentum k and the length scale 
j3 of the potential. The maximum value of k in the figures 
corresponds to an energy E — Uq. Left and right columns 
refer to different potential intensities: (a)-(c) a = 0.2, (b)-(d) 
a = 1. The plot in Figs. 18191 must be compared with the 
case /3 = 1 in (a) and (c) respectively. Dark regions indicate 
complete reflection or transmission, light gray (color online) 
corresponds to a 50% transparency. 



across the potential wells during the expansion. 

In the presence of blue-detuned speckles, see Fig. [5J 
the behavior is similar although in this case the conden- 
sate may undergo a reflection from the highest barriers 
that eventually stop the expansion as happens at the left 
side of the particular disorder realization in the figure. 
Even in this case the central part of the condensate gets 
localized, being trapped by two barriers that act as a po- 
tential well in the previous case Q . We have also verified 
that, as one would expect, the case of a gaussian random 
disorder is characterized by an intermediate behavior be- 
tween the former two, with part that is reflected by the 
highest barriers and part that is localized in the central 
wells. 

A central question is whether the observed behavior 
is of classical or quantum nature. Indeed, to observe 
non trivial localization effects caused by multiple inter- 
ference of the condensate in the speckle potential, the 
single wells/barriers should behave as quantum reflectors 
[25j . A qualitative insight on the behavior of the random 
potential can be therefore obtained by considering the 
case of a single defect. In case of the speckles a suitable 
model is a sech-squared potential of the form 



U (z) = aUosech 
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(12) 



where a and (3 are scaling factors for energies and lengths 
respectively. The transmission coefficient of this poten- 
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(13) 

with rj = 2mUoL 2 ) /Ti 2 and k = kL a = \/2mELo/Ji, E 
being the energy of the incoming wavepackct. For con- 
venience here we set the energy scale to the TF chemical 
potential of the condensate, Uq = htf- The length scale 
instead is fixed to Lq = ^l Cl with 7 ~ 2.72, in order to 
match, for = 1, the correlation length of the potential 
in Eq. 1|12[) with that of the random potential. With this 
choice the correspondence with the cases shown in Figs. 
15151 is for a = 0.2 and = 1. 

The ability of the above potential to act as a quantum 
reflector can be suitably quantified by introducing the 
function q(k,a,0) = 2 | 0.5 — T(k, a,0)\ — 1, that van- 
ishes in case of complete transmission or reflection, and 
equals one for a 50% transparency. In Fig. ^|we show a 
density plot of q as a function of k and for two values 
of a, considering both the case of a potential well and 
of a barrier (that are relevant respectively for the com- 
parison with the red and blue detuned speckles). The 
range chosen for the incident momentum k corresponds 
to energies up to Uq. The figure shows that in case of the 
speckles with a correlation length as in the experiment 
(0 = 1, a = 0.2) the range of energies where quantum ef- 
fects are evident is just a very narrow region close to the 
top of the barrier or at the well border. For this value of 
the correlation length even increasing the intensity of the 
potential (a factor of 5 in the figure) the situation does 
not change substantially. Instead, by reducing the length 
scale of the disorder (/3 — ► 0) quantum effects may even- 
tually become predominant in a wide range of energies. 
This corresponds to the fact that the height of the single 
defect should vary by a quantity at least of the order of 
the energy E of the incoming wave packet in a distance 
short compared to its de Broglie wavelength X,ib, that 
is: \dU/dz\ XdB > E. As discussed and experimentally 
demonstrated in 0, the above condition becomes very 
difficult to fulfil when the defects are created by near- 
infrared light as for a speckle potential. 

These considerations suggest the interpretation of the 
observed localization as a classical effect due to the actual 
shape of the potential. In this picture the condensategets 
partially localized by the presence of high barriers [5j or 
deep wells Q in the potential that act as single traps 
when the local chemical potential becomes of the order 
of their height. This is also confirmed by the comparison 
with the case of the periodic potential, that presents a 
qualitatively similar behavior as shown in Fig. 1111 Even 
in this case the most energetic part of the condensate 
expands nearly as free, whereas the bulk remains trapped 
in the central wells of the potential. The same picture 
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FIG. 11: (Color online) (a)-(d) Density profiles of the conden- 
sate (red continuous line) during the expansion in the waveg- 
uide in the presence of a periodic lattice (shown in frame (e)), 
for different times (t = 0, 25, 50, 75 ms, from (a) to (d)), com- 
pared with the free expansion case (blue dashed line). Notice 
that the y-sxis scale changes from (a) to (d). 



disorder have been also discussed. 

In the small amplitude regime the dipole and 
quadrupole modes are undamped and characterized by 
uncorrelated frequency shifts that depend on the partic- 
ular realization of disorder. This behavior, predicted by 
a perturbative sum rules approach, has been confirmed 
by the direct solution of the GP equation and observed 
in the experiment . The theoretical analysis shows also 
that the average features do not depend crucially on the 
particular kind of disorder, but are however significantly 
different from the periodic case. 

When released in a ID waveguide the condensate may 
be trapped into single wells or between barriers of the 
random potential, yielding a reduced expansion. These 
phenomena are of classical nature and take place prefer- 
ably near the trap center where the less energetic atoms 
reside. The outer part of the condensate instead expands 
almost freely, unless it encounters a high enough (reflect- 
ing) barrier. This behavior has been observed in recent 
experiments where the condensate is let expand in the 
presence of potential wells or barriers [f| . In the first 
case the qualitative behavior is very similar to that of a 
periodic system or even of a single well. 

We notice that in order to observe Anderson localiza- 
tion or related phenomena in a ID waveguide one should 
instead have interference of multiple quantum reflections 
of matterwaves. This regime could be achieved by reduc- 
ing the correlation length of the random potential, but 
may be not a trivial task due to the diffraction limit on 
the size of the defect created by light @|. 

In this respect the present analysis, besides provid- 
ing useful informations on the superfluid behavior of a 
condensate in the presence of a rough surface potential, 
suggests that it would be interesting to engineer other 
kinds of potentials by reducing the spacing or increasing 
the steepness. 



holds even in case of a single well, as discussed in |(|. 

Concerning the role of interactions, we note that they 
introduce the dephasing at the origin of the fast density 
modulations shown in the figures, that may eventually 
lead to a breakdown of the superfluid flow as discussed 
in Section Till Bl Their possible contribution to localiza- 
tion instead is not evident. Rather they act against local- 
ization, since they are responsible of the fast expansion 
of the lateral wings (the expansion in the noninteracting 
case would be much slower). 

V. DISCUSSION AND CONCLUSIONS 

A general analysis of the effects of a weak disorder cre- 
ated by speckle light on the collective modes and the 
expansion of an harmonically trapped condensate has 
been presented by using the Gross-Pitaevskii (GP) the- 
ory. The effects of different kinds of random potentials 
and a systematic comparison with the case of a periodic 
lattice with spacing of the order of the length scale of the 
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APPENDIX A: RANDOM DISTRIBUTIONS 

In this section we discuss how the random distributions 
used in the paper are constructed and characterized. For 
simplicity here we will use dimensionless units (lengths 
are expressed in units of an arbitrary scale £ whose actual 
value is irrelevant here). 

Following |27| the speckle distribution is constructed 
by starting from a random complex field <p(x) (on a grid) 
whose real and imaginary part are obtained from two 
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FIG. 12: (Color online) Continuous line (red): autocorrela- 
tion function (top) and intensity distribution (bottom) for the 
speckle potential in Fig. Q The (blue) dashed lines represent 
the expected average values over several realizations. 
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FIG. 13: (Color online) Autocorrelation function (top) and 
intensity distribution (bottom) for the gaussian random po- 
tential in Fig. Q The (blue) dashed lines represent the ex- 
pected average values over several realizations. 
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FIG. 14: (Color online) Autocorrelation function (top) and 
intensity distribution (bottom) for the periodic potential in 
Fig-0The (blue) dashed lines represent the expected average 
values over several realizations. 



independent gaussian random distribution rj(x) with zero 
mean (ri(x)} —0, unit standard deviation, and correlation 
function (ri(x)r)(y)) ~ 6(x — y). The speckle intensity 
field is then defined as 



where the operator J- indicates the Fourier transform 

F[<p\ = J dxtp(x)e 2mxy (A2) 
and W(y) indicates the aperture function 



W(y) 



1 i£\y\<D/2 
elsewhere 



(A3) 



The resulting distribution probability of the speckle in- 
tensities is 281 



P(I) = c 



- .-//(/) 



/</) 



(A4) 



and can be further normalized to 07 = (I) = 1 (the 
normalized speckle distribution is indicated in the text 
as vs(x)). The spatial (auto)correlation is 

T(Ax) = (I(x)I(x + Ax)} = 1 + sincpAx) 2 (A5) 

(the average (••) stands for an integration over x and 
an average over many realizations) with sinc(a;) = 
sin(7ra;)/(7ra;). The correlation properties can be sum- 
marized by the correlation length l c defined as the width 
at the half value of the maximum of T(Ax) (in Ax = 0) 
with respect to the background. In case of a one dimen- 
sional speckle distribution as that considered here l c is 
related to the aperture width by l c — 0.88/D. 

As a second source of disorder we consider a gaussian 
random distribution defined by |29| 



g(x)=T- 1 VTOllW 



(A6) 



I(x) = \T- X [W[y)T[ip{x) 



(Al) 



where r)(y) itself is a gaussian random distribution (de- 
fined as above) and the aperture function is W(x) — 
exp(— x 2 /2a 2 ). By using the properties of r] is then easy 
to demonstrate that both the real and imaginary part 
of g(x) are gaussian random distributions with a corre- 
lation function T(Ax) = W(Ax) and correlation length 
l c = 2\/21n2o-. Here we will consider in particular the 
imaginary component, I(x) = Im[g(x)] (in the text the 
normalized distribution is indicated as vg(x)). 

Finally let us discuss how to chose the wavevector of 
the periodic potential vl(x). In this case, the periodic- 
ity of the potential reflects in the periodic structure of 
the correlation function T(Ax) ~ 1 + 2cos 2 (gAx). By 
restricting over a single period the correlation length can 
be defined as above, and a straightforward calculation 
yields g — n/2l c . A suitable choice to compare the effect 
of disorder to the case of an ordered lattice described by 
the periodic potential vj,{x) is therefore to require the 
two potential to have the same correlation length. This 
seems a reasonable choice as shown in Fig. ^ 

The autocorrelation functions and the intensity distri- 
butions of the three potentials are shown in Figs. 1121141 
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